Coherence properties of the two-dimensional Bose-Einstein condensate 
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We present a detailed finite-temperature Hartree-Fock-Bogoliubov (HFB) treatment of the two- 
dimensional trapped Bose gas. We highlight the numerical methods required to obtain solutions 
to the HFB equations within the Popov approximation, the derivation of which we outline. This 
method has previously been applied successfully to the three-dimensional case and we focus on the 
unique features of the system which are due to its reduced dimensionality. These can be found in the 
spectrum of low-lying excitations and in the coherence properties. We calculate the Bragg response 
and the coherence length within the condensate in analogy with experiments performed in the quasi- 
one-dimensional regime [Richard et al., Phys. Rev. Lett. 91, 010405 (2003)] and compare to results 
calculated for the one-dimensional case. We then make predictions for the experimental observation 
of the quasicondensate phase via Bragg spectroscopy in the quasi-two-dimensional regime. 
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I. INTRODUCTION 

Bose-Einstein condensation (BEC) in (quasi-) two- 
dimensional systems has only recently been obtained in 
the laboratory [1, 2]. Thus, many properties have yet to 
be explored both experimentally and theoretically. We 
present an investigation of an isotropic two-dimensional 
BEC with the aim of providing detailed predictions for 
comparison with future experiments. 

The manner in which dimensionality can fundamen- 
tally alter the physics of a system is clearly apparent in 
the Mermin-Wagner-Hohenberg theorem, which forbids 
a spontaneously broken symmetry with long range or- 
der in a homogeneous two-dimensional system [3-5]. In 
terms of the coherence function G^ l '(x, x'), this means 
that lim\ x _ x i^ oc G^(x,x') ^ 0, which can be seen as 
the definition of BEC [6, 7], is impossible for T > 
in a uniform two-dimensional system. Thus, in a two- 
dimensional Bose gas BEC cannot occur at finite tem- 
peratures. Phase fluctuations make the formation of a 
globally coherent phase impossible. Despite this, a dif- 
ferent transition of the Kosterlitz-Thouless (KT) type [8- 
10] to a state with an analytical decay in the coherence 
function is possible in the ideal system. 

With confinement in a harmonic trap, the modified 
density of states allows the 2D system to Bose condense. 
Nevertheless, below the critical temperature there is a 
large phase fluctuating regime in which the superfluid is 
best described as a quasicondensate [11]. Unlike a true 
BEC, phase coherence only extends over regions of a size 
smaller than the extent of the condensate, characterized 
by the coherence length. This regime has been referred 
to as the KT phase, although the physical state of the 
interacting system in this phase fluctuating regime has 
yet to be thoroughly investigated. Phase fluctuations can 
enter the uniform gas in the form of vortex/antivortex 



pairs, or topological charges, which unbind at the point of 
the KT transition. Thus, the phase fluctuating state may 
well be a regular lattice of pairs of opposite topological 
charges in the sense of the KT phase, but this is not the 
only possibility and further investigation is required. 

In a previous publication [12] we have discussed how 
the semi-classical approximation fails to describe BEC 
consistently in two dimensions and have shown results to 
prove that these problems can be removed by applying 
the more complex Hartree-Fock-Bogoliubov (HFB) for- 
malism. The aim of the present publication is to present 
a detailed and more complete discussion of the properties 
of a two-dimensional BEC as is possible within the HFB- 
Popov approach. Our emphasis lies on the coherence 
properties which are crucial for the question of whether 
the superfluid state is best described as a BEC or as 
a quasicondensate. In the following section we outline 
the HFB formalism and explain our methods of obtain- 
ing solutions. Then, in Section III, we present our re- 
sults, such as the density profile of the condensate and 
non-condensate, the excitation spectrum and the coher- 
ence function. In Section III C we present the momen- 
tum profile and coherence length of a phase fluctuating 
condensate, indicating how these could be measured in 
forthcoming experiments. Our work is concluded in Sec- 
tion IV. 



II. FORMALISM 

A. Mean-field theory and HFB-Popov equations 

The time-independent, second quantized form of the 
grand-canonical many-body Hamilton operator for our 
system is given by 

H = jd 2 r ^(r) (h(r) - fi\ 4>(r) 
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(1) 



Here, h(r) = —j^A + f7 t rap( r ) is the single particle 
Hamiltonian with the external potential t/trap of the 
atom trap, and g is the coupling parameter that char- 
acterizes interparticle scattering. For collision processes, 
we assume a hard sphere potential within the usual 
pseudo-potential approximation [13], i.e. V(r — r') = 
(?<5( 2 )(r — r'). For a dilute gas this is a good approxi- 
mation, however care must be taken in determining the 
coupling constant g. Usually it is derived from an ap- 
proximation to the two-body T-matrix in the zero-energy 
and zero- momentum limit, as appropriate for scattering 
processes in an ultra-cold system. In three dimensions, 
the two-body T-matrix for a dilute gas is well described 
within the s-wave approximation, g — Airh a^D/m, where 
(23D is the s-wave scattering length. In two dimensions, 
however, the two-body T-matrix vanishes at zero energy 
[14]. Therefore, many-body effects introduced by the 
surrounding medium must be taken into account when 
studying two-dimensional gases. For a trapped gas, this 
leads to a spatially dependent coupling parameter g{v). 
Furthermore, the exact form of the coupling strength de- 
pends on the tightness of the confinement in the axial 
direction. With the parameters from [1], using the ter- 
minology of [14], we consider this system to be in the 
quasi-2D regime. Therefore, for the calculations under- 
taken in this work, we use the following approximation 
to the many-body T-matrix at zero temperature for the 
coupling parameter [14]: 
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The scattering length a2D in the quasi-2D regime is given 
by a 2B =AyJV/Bl z e~^ l *l a3D , B w 0.915. This result 
was first obtained by Petrov et al. [15, 16] by considering 
the 2D scattering problem. We will present a detailed 
study of interactions in the 2D Bose condensed system 
elsewhere [17]. 

We decompose the Bose field operators, in the stan- 
dard fashion [18, 19], into classical and fluctuation parts, 
?/>(r) ~ (V'(r)) + Stp(r) = ^o(r) + <5?/>(r), where the con- 
densate wave function ^o( r ) is normalized to the number 
of particles in the ground state, i. e. Jd 2 r |*o( r )| 2 = -^o- 
The Hamiltonian (1) can then be diagonalized by a uni- 
tary transformation to the quasiparticle operators on 
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d 2 r * (r) h(r) - /i + - g(r) n c (r) * (r) 
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(3) 



where C = h(r) — /i + 2g(r)n(r) and n c (r), n(r) and 
n(r) = n c (r) + n(r) are the condensate, non-condensate 
and total densities, respectively. The functions m, Vi are 
referred to as quasiparticle amplitudes, and Ei arc the 



quasiparticle energies. The first term in (3) is the con- 
densate part and merely a c-numbcr. The second term 
is the Hamiltonian for non-interacting quasiparticles and 
is formally equivalent to the case of the harmonic oscilla- 
tor. The constant energy shift C arises from the Bogoli- 
ubov transformation [18] and from terms left over from 
the quartet operator averages of the fluctuation operators 
which are factorized in a fashion analogous to Wick's the- 
orem [20, §4.2]. However, this energy shift has no impact 
on the solution of the HFB equations. 

The form (3) of the Hamiltonian requires that the order 
parameter obeys the generalized Gross-Pitacvskii equa- 
tion (GPE) 



h{v) - fA *o(r) + <?(r) (n c (r) + 2 n(r) ) * (r) - (4) 



and that the quasiparticle amplitudes and energies obey 
the coupled Bogoliubov-de Gennes (BdG) equations 



£u 4 (r)- 3 (r)*o(r) 2 w 4 (r) = E iUi (r) 
tvi{r) - g(r)%{rf Ui (r) = -E iVi (r) , 



(5) 



so as to eliminate off-diagonal terms in the quasiparti- 
cle field operators. The BdG equations determine the 
elementary excitation modes of the condensate. We re- 
fer to (4), together with (5), as the HFB equations. 
Note that we have taken the Popov approximation by 
neglecting the anomalous average of the fluctuation op- 
erator, m(r) = (5ip(r)5tp(r)}, whereby avoiding diver- 
gence problems of this quantity and the occurrence of 
a gap in the excitations spectrum [18, 19]. Once the 
BdG equations are solved, the non-condensate density 
h = (5^ (r)5ip(r)) can be obtained by populating the 
quasiparticle states, 
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where the quasiparticle distribution function with the in- 
verse temperature /3 is given by 



f B (Ei) = (a]a t ) 
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(7) 



Here, the fugacity z is determined by the difference be- 
tween the chemical potential /i and the condensate eigen- 
value A, z = e^^~^', since the quasiparticle energies are 
measured relative to the condensate [20]. To a good ap- 
proximation, we can use the result for the non-interacting 
gas, i. e. 
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The system we consider has a finite number of atoms. 
The fugacity fulfills the practical purpose of preventing 
the number of thermal atoms from exceeding the total 
atom number and, hence, the condensate density from 
becoming negative in our numerical calculations. 



In order to study coherence properties, we calculate the 
normalized first order correlation, or coherence function, 
which can be written in terms of the field operators as 

[21] 



3 (1) (r,r') 
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Given the decomposition of the field operator, the co- 
herence function can be expressed in terms of the off- 
diagonal densities 



n c (r,r') = *o(r)*o(r') 
n(r,r') = (^(r^r')) 



(10) 
(11) 



The latter can be calculated from the off-diagonal version 
of (6). Using the above for the correlation function, (9) 
gives 



(i), m _ n e (r,r')+n(r,r') 



>(r,r') = 



i/n(r) n(r') 
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The correlation function is related to the momentum 
spectrum of the condensate by a simple Fourier trans- 
formation, i. e. 



i(k) = (0t(k)^(k)> = / d 2 rdV e lk ' (r - r ' 



) /,/lt 



(^(r)V(r')) , 

(13) 

with 4>(k) and $ (k) being the field operators in momen- 
tum space. This implies that coherence properties can be 
directly measured in an experiment, as has been done in 
[22] for the quasi-one-dimcnsional case by means of Bragg 
spectroscopy. In a Bragg experiment, the propagation 
speed of the light field is determined by the detuning of 
the crossed laser beams [23]. Therefore, the spectral re- 
sponse of the condensate is measured as a function of the 
detuning. To establish the relationship with the momen- 
tum distribution, we use the relation between the detun- 
ing 8 and the momentum within the condensate plane p± 
for a n- photon process, 



5 



nk L p± 
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where k^ = 2ir/\, A is the wavelength of the light field 
(780.02 nm for Rubidium [22], 589 nm for Sodium [24]), 
and to the mass of the atoms. 



u±, Utmp = m{J\r 2 /2. Thus, our system is cylindri- 
cally symmetric and we can effectively treat the problem 
as one-dimensional upon changing to cylindrical coordi- 
nates. We scale all equations to computational units, 
i. e. lengths by the oscillator length, a,Q — y / Ti/muj±, and 
energies by the Rydberg of energy, Eo = Huj±/2. 

The calculation follows a self-consistent, iterative 
scheme, as proposed in [18]. First, the GPE is solved 
with the non-condensate density set to zero. Taking this 
calculated condensate density, the BdG equations are 
solved to obtain the quasiparticle modes. These are then 
populated through the quasiparticle distribution function 
(7), with the sum of all the excited particles yielding the 
thermal density. With the non-condensate density now 
known, we go back and solve the generalized GPE and 
the whole process is repeated until convergence. 

To begin, we expand the order parameter in a set of 
basis states. A convenient basis for this problem is given 
by the eigenstates of the 2D harmonic oscillator, since the 
single-particle Hamiltonian is diagonal in this basis. To 
take the cylindrical symmetry into account, we write the 
eigenfunctions of the oscillator problem h osc = — A + r 2 
in terms of the Laguerre polynomials L™, 
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r ™ e --+™v L™(r 2 ) 



(15) 

with the eigenenergies E n ^ m = 2(2n + to + 1). The quan- 
tum number to defines the angular momentum. Since the 
condensate ground state has zero angular momentum, for 
the solution of the GPE merely the to = subspace must 
be considered. In order to numerically solve the GPE, we 
use an optimization routine with a Thomas-Fermi profile 
as the initial guess. 

The solution of the BdG equations follows the method 
described in [19]. In a first step, the BdG equations (5) 
are decoupled by a transformation to the auxiliary func- 
tions ipf^ 0") = Ui (r) ± Vi (r) . Omitting spatial dependen- 
cies, this leads to 






M ) 2 Vf + 2gn c (h GP - M ) ^ +) = E 2 ^ 



v2 , 



/x) Vf ' + 2$(fc OP - mH ^V = Ef $ 



(16) 



where h GP = h + g (n c + 2n) is the Hamiltonian in the 
generalized GPE (4). The auxiliary functions V'i ±) ( r ) are 
then expanded in the basis set in which h GP is diagonal. 
This basis we term the Hartrce-Fock (HF) basis and it 
is obtained by the full solution of the generalized GPE, 



B. Numerical Methods 

We discuss some aspects important to the solution 
of the finite temperature HFB equations. The trap- 
ping frequency in the axial direction is sufficiently large 
so that the dynamics in this dimension are frozen out 
(hu> z > k^T). In the radial plane, we consider an 
isotropic trapping potential with the radial frequency 



(h GP - nj <f> a (r) = e a 4> a (r) , 



(17) 



where {4> a {r)} is the HF basis. The primary advantage 
of this basis is that all excited states are by definition 
orthogonal to the condensate, after the lowest momen- 
tum state, which is the condensate state itself, has been 
removed from the basis set. Both the calculation of the 



HF basis, as well as the solution of the decoupled BdG 
equations arc only linear problems, since the condensate 
density is given from the solution of the GPE, and can 
be solved in a straightforward manner. The eigenvalue 
problem corresponding to the BdG equations is block- 
diagonal with no overlap between the subspaces of dif- 
ferent angular momentum, so that the solution to (16) 
can be obtained separately in each subspace. The ther- 
mal density then follows from (6) by summing up the 
contributions from all angular momentum subspaces. 

Naturally, the number of basis states used in the dis- 
crete, quantum mechanical calculation is limited by an 
upper energy cutoff, e cu t, which must be introduced con- 
sistently in all angular momentum subspaces. To account 
for the contributions above the energy cutoff, we use the 
semi-classical approximation [25, 26], so that 

/>oo 

n{r) = ]T nf» * e(e cut - E t ) + / dE h sc (E, r) . 

(18) 

The contribution nf m (r) below the cutoff is given by the 
addent in (6), and above the cutoff by the semi-classical 
equation, with the Heaviside function 0, and again omit- 
ting spatial dependencies, 



15 ' * 

; \ 



10 



<= 







10 



15 



FIG. 1: Non-condensate density at T/T c = 0, 0.1, 0.25, 0.5 
and 0.75 (from bottom to top). The lowest line corresponds 
to the quantum depletion. 
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FIG. 2: Condensate (solid) and non-condensate density 
(dashed) at 0.5 and 0.9 T c . 
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III. RESULTS 



We will now present the results of our numerical calcu- 
lation. We consider a sample of 2000 sodium atoms that 
are trapped in a harmonic potential with the parame- 
ters of the experiment by Gorlitz et al. [1]. The radial 
trapping frequency is u± — 2ir x 790 Hz. Unless other- 
wise stated, all quantities are expressed in dimensionless 
form. 



A. Thermal density and condensate population 

1. Density profiles 

Figure 1 shows the thermal density at different tem- 
peratures. The temperature dependent term in (6) leads 
to the formation of the characteristic off-center peak of 
the non-condensate density. It is located at the edge of 
the condensate due to the repulsion of the thermal atoms 
by the condensate. This is depicted in Figure 2, where 
the two densities n c and n are plotted together. In com- 
parison to the rapidly decaying condensate density, the 
thermal density has a long tail. Thus, the condensate is 



relatively dense with a sharp peak within the diffuse ther- 
mal cloud. The tail of the thermal cloud becomes longer 
as the temperature increases, while the condensate radius 
does not change significantly even if the number of con- 
densate atoms drops by an order of magnitude. Note that 
the long tail contains a large number of atoms despite 
its low density because the spatial integral is weighted 
by a factor of r (r 2 in three dimensions). The slight 
change in the size of the condensate can also be seen 
in the shift of the non-condensate peak towards the trap 
center with increasing temperature, remaining located at 
the edge of the condensate. Note that even at zero tem- 
perature there is a small fraction of excited atoms due to 
the temperature-independent quantum depletion term in 
(6). 

From the existence of a well defined condensate and 
non-condensate density we can already infer coherence 
information. Following the argument in [12, 16], a mea- 
sure of phase fluctuations is given by (S 2 ) « h/n c , where 
5 is the phase fluctuation operator in the alternative de- 
composition tp( r ) — \/n(r) e l5 ( r \ Phase fluctuations be- 
come important when (S 2 ) > 1. Thus, as long as n c > h, 
these fluctuations are suppressed in the system. 



2. Ground state population 

In Figure 3 the condensate population is shown as a 
function of temperature. The results are compared to 
the case of the trapped ideal gas, where the population 



is determined by a power law expression. We fit the 
following functional form to the numerical data: 



No =1 _(T_ 



where T r = aT r 



(20) 



In the case of the ideal gas, (3 — 2 and a = 1. In the fit 
the critical temperature is reduced by a factor of about 
5% with a ~ 0.95. This shift has two contributions: The 
finite size of the system reduces the critical temperature 
[27], but it is also modified by the interactions. This 
second contribution as been extensively discussed in the 
literature, see the recent publication [28] and references 
therein. For the exponent we find (3 « 1.70, which is 
15% smaller than for the ideal trapped gas. With these 
values, (20) parameterizes our data very well except near 
the critical temperature, where finite size effects are sig- 
nificant and the exact method of determining the shift 
of the chemical potential from the condensate eigenvalue 
becomes important. In the ideal gas the chemical po- 
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FIG. 3: Condensate population versus temperature. The solid 
line corresponds to a fit to (20), the dashed line shows the 
ideal gas power law dependence. The points are results from 
the HFB-Popov calculation. 

tential is zero at the point of the phase transition, and 
therefore the transition point is strictly defined. In the 
interacting gas, the chemical potential depends implic- 
itly on the non-condensate [19] and, when this becomes 
large, the transition point becomes smeared out. Techni- 
cally, the occurrence of the finite temperature tail can be 
explained through the fugacity term in the quasiparticle 
distribution function (7), which is explicitly given by (8). 
The term oc 1/Nq prevents the condensate population 
from becoming negative. However, this expression for 
the fugacity is only approximate and the shape of the tail 
and the speed with which it approaches zero depends on 
the explicit choice of the fugacity term at temperatures 
around T r . 



B. Condensate excitations 

The low-lying collective or elementary excitation 
modes of the condensate, determined by the solution of 
the BdG equations, are of interest because they reflect 
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FIG. 4: Low-lying excitation modes of the condensate as a 
function of temperature. The uppermost line corresponds to 
the breathing mode with angular momentum quantum number 
m — 0, the middle line to the quadrupole mode with m = 
2. The lowest line is the Kohn mode, m — 1, which lies 
constantly at the trapping frequency. 



certain fundamental symmetry properties of the system, 
as well as being easily accessible to experiments. 

Figure 4 shows the modes with angular momentum 
to = 0,1,2 as a function of temperature. Each of the 
three branches corresponds to the lowest quasiparticle 
energy eigenvalue in the lowest three, separated, angular 
momentum subspaces, in which the BdG equations are 
solved, c. f. Section II B. 

a. Breathing mode. The breathing mode corre- 
sponds to an oscillation of the condensate radius and 
lies at a frequency that is twice the trapping frequency. 
As shown in [29], this is due to a hidden symmetry of 
the many-body Hamiltonian with a <5( 2 )-interaction po- 
tential and a harmonic trapping potential in two dimen- 
sions. As the temperature increases, the non-condensate 
density grows and starts to constitute a deviation from 
the harmonic oscillator potential in the effective potential 
of the Gross-Pitaevskh equation, so that the frequency 
shifts slightly from 2ui±. The effective potential is weak- 
ened by the presence of the static thermal atoms so that 
the frequency decreases. If the dynamics of the thermal 
cloud were included in the calculation [30], then the full 
symmetry of the Hamiltonian would be restored and the 
mode frequency would remain precisely at 2Hlo±. 

b. Kohn mode. The Kohn mode corresponds to a 
center of mass oscillation of the whole condensate. Less 
effected by the perturbation to the harmonic potential 
of the static thermal cloud, it remains very constant at 
the trapping frequency, as is predicted by the generalized 
Kohn theorem [31]. However, looking closer at Figure 4, 
one may see a slight increase in the frequency near the 
critical temperature as the effective potential becomes 
less harmonic and, therefore, breaks the Kohn theorem. 
In our calculation we treat the thermal cloud as station- 
ary. An inclusion of the full dynamics of the thermal 
cloud would, again, ensure the Kohn mode remains con- 
stant at all temperatures [19, 30, 32]. 



c. Quadrupole mode. The quadrupole mode is the 
only low-lying mode which depends strongly upon the 
temperature. The frequency of this mode could thus, in 
principle, be used as a measure of the temperature of the 
2D gas. 

With increasing temperature, all three frequencies 
smoothly approach the frequencies of the non- interacting 
gas and the breathing and quadrupole mode become de- 
generate. Using a local density approximation with the 
relation jjl = g[n c (r) + 2n(r)] for the uniform gas in the 
Hartree-Fock approximation [33, §8.3], it is easy to show 
that the BdG equations for n c (r) = 0, given by 

(h(r)-fi + 2gh(r) - E^ u t (r) = , (21) 

recover the energies of the harmonic oscillator. In the 
case that there is still a condensate, the total density in 
the region where n c (r) =/= is approximately constant 
just below the critical temperature, so that the mean- 
field energy only constitutes a near constant shift to the 
trapping potential and, hence, only slightly alters the 
eigenfrequencies of the trap. 

We would briefly like to draw comparison with the 
three-dimensional case where the frequency spectrum 
looks very similar [19, 34]. The striking difference is the 
breathing mode which is temperature dependent in three 
dimensions, whereas it is a feature of the two-dimensional 
system to have breathing oscillations with a universal en- 
ergy of 2hoj± . 
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FIG. 5: Correlation function g( ' for the non-interacting Bose 
gas at different temperatures T/T c : 0, 0.025, 0.05, 0.1, 0.5, 
0.9 and 1 from right to left. 



condensate that smoothes the sharp corner as g^ drops 
to zero, causing a loss of numerical accuracy as we divide 
two very similar small numbers in (12). 

The coherence length can be extracted by measuring 
the full width at half maximum (fwhm) of g 1 - 1 ' and is 
shown in Figure 8. 

b. Non-interacting gas. The off-diagonal density 
matrix is known in closed analytical form for the non- 
interacting gas. In can be determined by means of the 
inverse Laplace transform of the zero-temperature Bloch 
density matrix [35, 36]. Its explicit form in two dimen- 
sions at temperature T is given by 



C. Coherence properties 

1. Correlation function 

a. Interacting gas. In Figure 5 the correlation func- 
tion gW(0,r) is depicted at various temperatures. The 
r-axis has been scaled by the size of the condensate. This 
is not the Thomas- Fermi radius, but we choose a minimal 
allowed condensate density in such a way that the whole 
condensate at zero temperature is phase coherent. 

The decay of the correlation function allows for a char- 
acterization of the gaseous system. At low temperatures 
the correlation function has a constant value throughout 
the extent of the condensate, indicating a truly coherent 
Bose-Einstein condensed phase with off-diagonal long- 
range order. Algebraic decay is associated with the KT 
phase and, at intermediate temperatures, the superfluid 
must be identified as a quasicondensate. At very high 
temperatures, clearly visible for the highest temperature 
in Figure 5, the coherence function decays exponentially, 
showing that long-range order is lost completely. 

At very low temperatures the correlation functions 
show some unphysical oscillations that are purely nu- 
merical noise. At low temperatures n(r) ps n c (r). In 
the limit h = the correlation function (9) is given by 
the Heaviside function 0(1 — r/r con ). However, there is 
a small contribution from the quantum depiction of the 
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We find the chemical potential /i = /i(T) for the 
trapped non- interacting gas by solving J2 n=a fB(E n — 
n + 1,/U, T) (n + 1) — TV = with respect to /i. Here, 

/b = [(l + JV - 1 )e^ £ "-' i )-l]" 1 is the Bose-Einstein 
distribution function with the fugacity factor (8) that 
takes into account the number of condensate particles 
from the HFB calculation. Without this fugacity factor, 
finite size effects, taken into account in the HFB calcula- 
tion, would be neglected and, therefore, the comparison 
would be between two approaches based on different as- 
sumptions. The expression (22) for the correlation func- 
tion is exact at all temperatures. Numerically, the infinite 
sum can be calculated up to any required accuracy. 

We compared our code against these exact results for 
the non-interacting gas. At all temperatures the HFB 
results agree perfectly with the correlation function cal- 
culated from (22), implying that the numerics works well 
even at high temperatures. Deviations would indicate 
an insufficiency in the basis set or inaccuracy due to an 
insufficient fineness or range of the computational grid. 

The influence of the interactions can be seen in Figure 




FIG. 6: Correlation function g^ 1 ' for the non-interacting Bose 
gas at different temperatures T/T c : 0.025, 0.1, 0.5, 0.8 and 
0.95 from right to left. The solid lines correspond to the 
HFB-Popov results for the interacting gas, the dotted lines 
represent the exact result for the non-interacting gas (22). 




FIG. 7: Momentum spectrum, with intensity corresponding 
to the fraction of scattered atoms as function of the Bragg 
detuning 8. Temperatures T/T c from top: 0, 0.025, 0.05, 0.1, 
0.5, 0.9 and 1. 



6. The solid lines show the coherence function for the 
interacting gas, compared to the exact non-interacting 
gas equation, shown as dotted lines. We see that interac- 
tions increase the coherence length in a large part of the 
temperature regime. At 0.8 T c and above, however, the 
coherence length of the interacting gas is decreased rela- 
tive to the non-interacting gas. This can be explained by 
the effect of the mean-field interaction on the condensate 
radius. At high temperatures (7V C small) the radius of 
the interacting condensate is approximately the same as 
the radius of the non-interacting condensate, given by the 
size of the lowest harmonic oscillator state. At low tem- 
peratures, however, the condensate population is large 
and mean-field effects broaden the condensate. Corre- 
spondingly, the coherence function of the interacting con- 
densate is broader. If the spatial coordinate was scaled 
by the size of the condensate as in Figure 5, the plot 
would show that interactions always reduce the range of 
coherence. 



2. Momentum spectrum 

Figure 7 shows the momentum spectrum correspond- 
ing to (13), as it could be measured by means of Bragg 
spectroscopy. It has been calculated by Fourier trans- 
forming the correlation function shown in Figure 5. On 
the ordinate is the detuning of the Bragg laser beams, 
which is directly proportional to the momentum of the 
atoms, c. f. (14). The highest peak corresponds to the 
lowest temperature where the momentum distribution of 
the atoms is narrowest. With increasing temperature the 
spectrum is broadened. An experimental setup is limited 
by its resolution at low temperatures, because of the de- 
crease of the spectral width. 



In Figure 8, the coherence length obtained from the 
momentum spectra in Figure 7 is plotted. We deter- 
mine the coherence length by fitting the momentum pro- 
file and measuring the half width of the fit. A Gaussian 
provides a good fit at temperatures <0.9T C . Above a 
small crossover regime, the momentum profiles at tem- 
peratures higher than 0.925 T c fit more closely to a Lo- 
tcntzian. Thus, the data points in Figure 8 correspond 
to the fwhm of a Gaussian or a Lorentzian fit, depending 
on which gives better agreement with the data. 

In the same graph the lengths are compared to those 
obtained from the results shown in Figure 5. Qualita- 
tively the results agree with each other, although those 
obtained from g 1 - 1 ' lead to somewhat smaller values for 
the coherence length. Also one can see that the extracted 
half widths of the correlation function are subject to a 
slight inaccuracy, whereas the widths calculated from the 
fits to the momentum profile result in a smooth line over 
the whole temperature regime. Note that the coherence 
length is to some extent a matter of personal definition, 
as e. g. we could have chosen 1/e rather than the fwhm. 

Looking at the decay of the coherence length in Figure 
8, we can distinguish three different regimes. Close to 
zero temperature, the slope is very steep and the coher- 
ence length decreases to a third of the condensate size 
by about 0.1 T c . Then, up to about 0.8 Tc, 1$ decreases 
monotonically, but much more slowly. From there up 
to the critical temperature, the coherence length again 
drops rapidly. 

A decreasing coherence length directly implies a loss 
in the global phase coherence of the superfluid phase. 
A true Bose-Einstein condensate cannot be said to exist 
when the phase of the order parameter fluctuates on a 
length scale significantly smaller than the extent of the 
condensate. At this point we should instead refer to a 
quasicondensate. From Figure 5 we see that the coher- 
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FIG. 8: Coherence length of the condensate, shown in both 
real units and scaled by the extension of the condensate. The 
upper line (+) has been calculated from the momentum spec- 
trum (left panel), the lower line (■) is the data from Figure 
5. 



ence length drops smoothly. Therefore, it is difficult to 
determine an exact point on the temperature scale where 
the transition from a true condensate to a quasicondcn- 
satc takes place. At about 0.5 T c the coherence length has 
dropped to about half of the maximal value. The max- 
imal value we can use to determine the spatial extent 
of the condensate, indicated on the right axis in Figure 
8. The treatment in [16] predicts a value of approxi- 
mately 0.4 T c for our parameters for phase fluctuations 
to become dominant. Our result is, therefore, consistent 
with [16], although the coherent phase seems to persist 
at slightly higher temperatures. 



3. Comparison to ID 

Similar behaviour has been observed in calculations 
for the one-dimensional Bose gas at finite temperatures 
[37]. Looking at the coherence function presented by 
Ghosh, we see that in the one-dimensional case the co- 
herence length drops even more rapidly than in the two- 
dimensional case, showing that phase fluctuations be- 
come much more dominant as the dimension is reduced 
further. In ID, the temperature range between 0.3 and 
0.5 T c has about the same coherence properties as the 
range around 0.9 T c in our 2D calculation. Ghosh iden- 
tifies the ID phase at temperatures as low as 0.1 T c as a 
quasicondensate with large phase fluctuations. From an 
examination of Figure 5, we see that, at this temperature 
in the 2D case, even if the coherence length has decreased 
slightly, there is still a large proportion of the condensate 
where g^ 1 ' is constantly 1, indicating that the system is 
essentially a phase coherent BEC. 

In ID the Lorcntzian momentum profile has been 
found to be characteristic of the phase-fluctuating 
quasicondensate [38] and has been used as an the 
identifying signature of such a phase [22]. However, 



we are convinced that the shape change we observe 
is not a signature of a phase fluctuating condensate, 
but the effect of the fugacity term as Nq goes to zero. 
Furthermore, from looking at the correlation function 
in Figure 5, we would expect the phase fluctuations to 
become important, indicating the presence of a quasi- 
condensate, at about 0.5 T c , long before the momentum 
profile becomes Lorentzian in character. 



IV. CONCLUSION 

We have used the HFB formalism to investigate the 
finite-temperature physics of a Bose-Einstein condensate 
confined to a two-dimensional geometry. Unlike the three 
dimensional case, phase fluctuations must be taken into 
consideration at comparatively lower temperatures. In a 
regime below the critical temperature they destroy the 
global coherence of the condensate and the superfluid 
state is best described as a quasicondensate. In the HFB 
formalism phase fluctuations arc included via the contri- 
bution to the non-condensate density from low-energy 
quasiparticles. We have shown that the formalism is 
not only applicable in the strictly phase coherent regime, 
but also that the quantities obtained, such as the single- 
particle off-diagonal density matrix, allow for a quantita- 
tive analysis even in the phase fluctuating regime. Our 
work is consistent with [16], although we find that, within 
the HFB treatment, the pure condensate phase persists 
to higher temperatures. 

The coherence length of the condensate can be deter- 
mined from its correlation function or the momentum 
profile. Following Aspect et al. for the one-dimensional 
case [22, 39], we have calculated the Bragg spectrum for 
a condensate in two dimensions. We found the values 
extracted for the coherence length to be in qualitative 
agreement with those calculated for the one-dimensional 
Bose gas, although a true BEC with global phase coher- 
ence still exists at much higher temperatures than in the 
ID case. The Bragg spectrum provides a clear signature 
of the quasicondensate phase and we anticipate experi- 
mental efforts in this area in the near future. 
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